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ABSTRACT 

After 25 years of quiescence, the microquasar V404 Cyg entered a new period of activity in June 2015. This X-ray source is known 
to undergo extremely bright and variable outbursts seen at all wavelengths. It is therefore an object of prime interest to understand the 
accretion-ejection connections. These can, however, only be probed through simultaneous observations at several wavelengths. We 
made use of the INTEGRAL instruments to obtain long, almost uninterrupted observations from 2015 June 20*, 15;50 UTC to June 
25*, 4:05 UTC, from the optical V-band, up to the soft y-rays. V404 Cyg was extremely variable in all bands, with the detection of 
18 flares with fluxes exceeding 6 Crab (20-40 keV) within 3 days. The flare recurrence can be as short as ~20 min from peak to peak. 
A model-independent analysis shows that the >6 Crab flares have a hard spectrum. A simple 10-400 keV spectral analysis of the 
off-flare and flare periods shows that the variation in intensity is likely to be due to variations of a cut-off power law component only. 
The optical flares seem to be at least of two different types: one occurring in simultaneity with the X-ray flares, the other showing a 
delay greater than 10 min. The former could be associated with X-ray reprocessing by either an accretion disk or the companion star. 
We suggest that the latter are associated with plasma ejections that have also been seen in radio. 
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1. Introduction 


V404 Cyg (hereafter V404) is a low mass X-ray binary (LMXB) 
consisting of a black hole (BH) with mass estimates ranging 
from ~9 to 15 Mq, and a O.T/g^Mo K3 III companion (Casares 
& Charles 1994; Shahbaz et al. 1994; Khargharia et al. 2010), 
located at a parallax distance 2.39+0.14kpc (Miller-Jones et al. 
2009). The inclination of the binary’s rotational axis is 67°*j 
(Shahbaz et al. 1994; Khargharia et al. 2010), the orbital period 
6.5 d (Casares et al. 1992). This transient underwent three peri¬ 
ods of outbursts during the 20th century (Richter 1989), the last, 
in May 1989, leading to its discovery as an X-ray transient by 
the Ginga satellite (as GS 2023+338, Makino et al. 1989). V404 
showed bright X-ray flares on short time scales (e.g. Makino 
et al. 1989; Terada et al. 1994), and this makes it an excellent 
source to study the connections between the accretion and ejec¬ 
tion phenomena, which are the probable origin of this behavior. 
V404 is one of the closest stellar-mass BHs, making it a rare 
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case where quiescence could be studied in detail, and variable 
remnant activity, attributed to a compact jet, was detected from 
radio to hard X-rays (e.g. Hynes et al. 2004; Xie et al. 2014). 
V404 is one of the few sources that defines the radio/X-ray cor¬ 
relation over a large range of luminosities, down into quiescence 
(Corbel et al. 2008). The good knowledge of the quiescent state 
makes understanding new outburst observations paramount as 
they allow the mechanisms responsible for the increased activity 
to be probed. 

On 2015 June 15 (MJD 57188), V404 went into outburst 
again. It was first detected by Swift (BAT and XRT) (Barthelmy 
et al. 2015), and then with MAXI, and INTEGRAL (Negoro et al. 
2015; Kuulkers et al. 2015). These early alerts triggered follow¬ 
up observations at all wavelengths. Preliminary results all report 
the detection of the source, variations of specific spectral fea¬ 
tures, and an extreme flaring activity at all wavelengths (e.g. 
Mooley et al. 2015; Munoz-Darias et al. 2015; Ferrigno et al. 
2015; Motta et al. 2015; Tetarenko et al. 2015a,b; Tsubono et al. 
2015). We triggered our INTEGRAL ToO program to obtain 
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Fig. 1. Left: INTEGRAL LCs of V404 in four spectral domains (a larger version of the plot including all energy ranges is available as online 
material). Right: 20-40 keV count rate (top), and 3-13/20-40 keV softness ratio (SR, bottom). The dashed horizontal line corresponds to 1 Crab 
(20-40 keV). MJD 57193 is 2015 June 20*^ 


quasi-continuous coverage in X/y-rays and in the optical V band. 
The first detection of multiple flares exceeding 30 Crab in 20- 
40keV, and possible correlated flaring in the V-band were re¬ 
ported in Rodriguez et al. (2015a), and Domingo et al. (2015). 

Our observations caught the source during the most intense 
and variable phase of this new outburst. Here, we hrst give the 
details of the observations and data reduction (Sec. 2). We then 
focus on the flaring behavior at high energies (Sec. 3) that we 
compare to the activity in the optical V-band (Sec. 4). A pre¬ 
liminary phenomenological spectral characterization of different 
intensity intervals is presented in Sec. 5. We discuss our results 
and compare V404 to other microquasars in Sec. 6. 


2. Observations and Data Reduction 

Our ToO program (Fig. 1) covered MJD 57193.66-57198.17 
(2015-06-20 15:50 UTC to 2015-06-25 4:05 UTC), i.e., IN¬ 
TEGRAL revolutions #1555 (continuous coverage) and #1556 
(two periods). The data of all the INTEGRAL instruments 
(see Winkler et al. 2003, and references therein for all instru¬ 
mental details) were reduced with the Off line Scientific 
Analysis (OSA) vlO.l software suite, with the latest calibra¬ 
tion files available at the time of writing. 

Images and 100 s binned light curves (LC) from the Joint 
European X-ray Monitors (JEM-X) and the Imager on Board the 
INTEGRAL Satellite (IBIS) were produced in two bands (3-13, 
and 13-30 keV) for JEM-X unit 1, and in four bands (20^0, 
40-80, 80-150, and 150-300keV) for the IBIS Soft Gamma- 
Ray Imager (ISGRI). 

The event data of the Spectrometer on INTEGRAL (SPI) have 
been fitted with models for the celestial sources and instrumental 
background following standard reduction processes. The 20-100 
keV LC of V404 as well as the other sources in the held were 
obtained in bins of 400 s. Background models have been built 
based on the pre-haring data of a representative empty sky re¬ 
gion, adjusting the normalization coefficient per hour (see, e.g., 
Strong et al. 2005, for a more general description of the method). 


As source intensity and hardness vary strongly on short 
time scales, we extracted luminosity/hardness dependent JEM- 
X, ISGRI, and SPI spectra over specihc time intervals of clean 
data. The spectra from the same time intervals were jointly ht- 
ted within XSPEC vl2.8.2. Since the instruments’ responses 
are possibly different for the high intensities observed, only 
phenomenological spectral hts are presented and the ht results 
should be viewed with some caution. 

The INTEGRALjopixcdl Monitoring Camera (OMC) huxes 
and magnitudes were derived from a photometric aperture of 
3x3 pixels (1 pix. =17.504"), slightly circularized, i.e. remov¬ 
ing one quarter pixel from each corner (standard output from 
OSA). The photometric aperture was centered on the source co¬ 
ordinates (default centroid algorithm) and did not include any 
signihcant contribution from other objects. We removed mea¬ 
surements with a severe problem flag, and, to restrict the noise, 
only measurements of 50 and 200 s duration were considered. 

3. Model independent description of the flaring 

Multi-wavelength LCs of V404 from the V band up to y-rays are 
highly structured with several large flares separated by calmer 
periods seen in all bands (Pig. 1, and see also Pig. 4 for a plot 
with all energy ranges). In the following, count rates (CR) are 
given in the ISGRI 20^0keV range. When the source CR in¬ 
creased above ~150-200 cts/s an intense X-ray flare systemati¬ 
cally followed. In the following, we thus set 1 Crab^ as the typ¬ 
ical limit between the off-flare and flaring intervals. We identify 
18 main events, i.e. peaks that reached at least 6 Crab (labeled 
with Roman numerals^ in Pig. 1, their main characteristics are 
given in the online table 1), with 11 exceeding 20 Crab during 

' The ISGRI/20^0 keV CR of the Crab is 165 cts/s 
F 2 o_ 40 kev =8 erg/cm^/s for a power law spectrum with r=2. 1 and 
a normalization of 10 ph/cm“^/s at 1 keV. 

^ V and XII contain two distinct events hardly distinguishable in Fig. 1. 
They appear under the same label in Fig. 1 (to keep it clear), and are 
named with a/b sub-labels in the text and table 1 . 
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Fig. 2. V-band (green), 3-13 keV (red), and 20-40 keV (blue) LCs 
around flares I (a) and X (b). The inserts show the cross-correlation 
functions of the 3-13 keV (red) and 20-40 keV (blue) vs the optical 
LCs over the same time intervals as the LCs. The dashed vertical lines 
represent the 0 lag level. 

our observations. Flares IV, XI, and XIII are the brightest we ob¬ 
served, reaching 43 Crab. The flares occur isolated (e.g. Ill, IV, 
VI) as well as in groups with peak to peak intervals as short as 
22 min (Va, Vb) The flares last 0.4-2.4 hr except peaks IV and 
XIII. The former shows a rather broad profile and has multiple 
peaks. This event lasted 4.8 hr in total and it is the longest flare 
of our observation. The latter reached about 40 Crab.The peak it¬ 
self lasted about 1.5 hr, but was preceded by a ~3 hr long, 3 Crab 
plateau seen only above 13 keV. It was followed by flares XIV 
and XV which show decreasing peak values. 

The 3-13/20^0 keV softness ratio (SR, Fig. 1, right) shows 
that the large variability of the source is associated with vari¬ 
ations of SR from ~0.03 to ~1.3 corresponding to F=0.1-2.5 
in simulated JEM-X/ISGRI power law spectra. SR~0.6 corre¬ 
sponds to F=2.0. Large spectral variations are visible in the off- 
flare intervals (Fig. 1, right). All the flares are hard, and they 
even all have SR<0.4 (oF<1.8). 


4. Optical vs. X-ray behavior 

The comparison of the optical (OMC) and X-ray (JEM-Xl and 
ISGRI) LCs shows a non trivial relationship. Significant flaring 
activity is evident in the V band EC (Eig. 1, left), with at least 
12 clear flares. The optical flare typically last 0.24-2.5 hr. While 
some events occur in simultaneity with X-ray flares, the optical 
emission is delayed with respect to the X-rays in other cases. 
Eig. 2 shows typical examples of these different behaviors. The 
cross-correlation function (ccf) between the X-ray and optical 
emission confirms the absence of lags for some of the flares (e.g. 
flare I that causes the peak at 0 in the ccf of Eig. 2a), and delayed 
optical emission from 1.5 min to 20-30 min is seen in others. 
The ccf of Eig. 2b shows an example of a ~3 min lag, while the 
ccf of Eig. 2a (in addition to showing the simultaneity of peak 
I and its optical counterpart) shows lag at 20-30 min represent¬ 
ing the delay between the small X-ray flares preceding peak I 
(around MJD 53193.9) and the subsequent optical flare (around 
MJD 53193.92). Due to the time resolution of OMC, however, 
lags shorter than 1 min can not be measured with our data, and 
additional lags of the order of seconds are not excluded. 

While most of the flares show a fast rise similar to the flares 


observed in X-rays, the two optical events occurring close to X- 
ray peaks IV and XIII seem to be exceptional. All other flares 
show fast rises (<1 hr), but these two events have slower rises 
(about 10 and 4 hr, respectively), and are both coincident with 
hard plateaus that precede the X-ray peaks. 

5. Spectral analysis 

We accumulated spectra from the brightest flares 
(CR>1000 cts/s) and the off-flare intervals (CR<165 cts/s). In 
the latter case we also only retained the hard intervals (SR<0.6) 
in order to exclude the softening visible after MJD 57197.9 
(Eig. 1). The resultant “v-fy” spectra are plotted in Eig. 3. 

The off-flare spectrum is well fitted (;^fy=1.2, 66 degrees 
of freedom, dof) by a model consisting of a power law with 
a high-energy cut-off dominating at 10-100 keV plus an 
additional power law dominating above 100 keV. The former 
component has F=1.0;|;g4, Ecut^lb;!;^ Efoid=23+5 keV, 

the latter has F=1.9 ^q 3, (F is the photon index defined as 
N(E)ocE“^.) Normalization constants were included to account 
for potential cross-calibration issues, or differences in the 
effective exposures (deadtime corrections, telemetry drop out). 
When set to 1 for ISGRI, we get ~1.9 for SPI and ~0.6 for 
JEM-Xl. The 20^00 keV (ISGRI) flux is -10^* erg/cm^/s, 
and the above model leads to an extrapolated 0.1-10^ keV flux 
~3.8 10“^ erg/cm^s, i.e. about 2% l^Edd for a 9 Mq BH. 

The flare spectrum is well represented (xl~0.9, 78 dof) by a 
single cut-off power law with F= 1.54+0.06, Ecut= 14.0!|^ keV, 
Efoid=87^j keV. An extra power law component is not sta¬ 
tistically required according to an E-test. The normalization 
constants are both close to 1.1 The 20^00 keV (ISGRI) flux 
is ~10“’erg/cm^/s, which leads to an extrapolated 0.1-10^ keV 
flux ~3 10“’erg/cm^/s, i.e. about 20% LEdd for a 9 M© BH. 

6. Discussion 

Over the 4 days covered by our INTEGRAL ToO, V404 showed 
a high level of emission with sporadic flares with a maximum 
20^0 keV dynamical range of 940 (flare XVI). During its flares 
V404 became the brightest X-ray object in the sky. In the hard 
off-flare state the spectral analysis shows the presence of two 
spectral components: a cut-off power law usually attributed to 
thermal Comptonization and an extra power law at energies be¬ 
yond 100 keV. Hard tails have now been seen in a large num¬ 
ber of systems (e.g. GRS 1915+105, Swift J1753.5-0127, GX 
339-4, or Cyg X-1, Rodriguez et al. 2008b; Tomsick et al. 2015; 
Joinet et al. 2007; Rodriguez et al. 2015b) and their origin is 
still highly debated, although a compact jet origin is favored in 
the case of Cyg X-1 (e.g. Russell & Shahbaz 2014; Rodriguez 
et al. 2015b). It is interesting to remark that the flaring activity 
seems primarily due to spectral variations of the cut-off power 
law only. We estimate an integrated luminosity L~ 0.2LEdd for 
the >6 Crab flares. This is a lower limit to the maximum lumi¬ 
nosity reached at the peak of the brightest flares. Eirst because 
we averaged the data, without isolating the brightest portions of 
the flares. Also we did not consider all the contributions below 
10 keV (disk, jet) that can provide a significant fraction of the 
bolometric luminosity. Assuming a simple scaling between the 
CR and LEdd with a constant shape for the variable component 
between the off-flare hard state and the flares, we conclude that 
all peaks with CR>3000 cts/s (~18 Crab) (flares III, IV, Va,b, 
VIII, X, XI, XIIa,b, XIII and XIV, Eig.l) reached -LEdd- 

The optical activity is also highly variable and shows flares 
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Fig. 3. y-fv spectra obtained from the fits to the off-fiare and flare inter¬ 
vals. The lines represent the spectral components used for the fits. 


(Fig. 1, 2). Some optical flares occur in conjunction with the X- 
rays, while other activity periods show delays (Fig. 2). The first 
specific length scale of this system is the separation between the 
BH and the companion: using the system parameters from Sec. 1 
and a 9 Mq BH we estimate ~2.2 lO'^ cm or ~75 light seconds. 
Hence, when no delay between the optical and the X-rays is ob¬ 
served (e.g flare I in Fig. 2a), the mechanism producing the op¬ 
tical emission could be related to X-ray reprocessing, either by 
an accretion disk or by the companion. The maximum delays 
expected would be around 60 s (outer disk), and <150 s (com¬ 
panion located at superior conjunction). 

Optical lags >10 min could be related to variable jet proper¬ 
ties, either as their intrinsic synchrotron emission, or from their 
interaction with the surrounding medium. Radio, and millime¬ 
ter (mm) flaring activity ascribed to discrete ejections have been 
reported during this outburst (e.g. Sivakoff et al. 2015; Mooley 
et al. 2015; Tetarenko et al. 2015b). Delays between X-rays and 
longer wavelengths are expected in the case of adiabatically ex¬ 
panding ejecta (van der Laan 1966; Mirabel et al. 1998), and we 
remark that Tetarenko et al. (2015b) report a mm peak on MJD 
57195.548, about 26 min after the MJD 57195.53 X-ray event (a 
<6 Crab flare that occurred slightly after flare XII). Interestingly, 
this delay also corresponds to the timescale of the mm flare in¬ 
crease (Tetarenko et al. 2015b) which renders the mm compati¬ 
ble with being causally related to the X-ray flare. As in other well 
known sources, e.g. XTE J1550-564, GRO J1655-40, GRS 
1915-1-105 (Fender 2006), in V404 discrete ejections may be 
causally connected to the X-ray activity. 

V404 reached maximum absolute peak intensities that are 
rather usual during outbursts of microquasars. If V404 was at 
10 kpc the maximum peak intensity would have been 2.5 Crab, 
a value similar to the aforementioned other microquasars during 
their brightest outburst(s) (e.g. Remillard & McClintock 2006). 
Short recurrent (multi-wavelength) flares have, however, been 
seen only in GRS 1915-1-105 (e.g. Greiner et al. 1996), and in 
V404 the flaring activity similarly recurs on timescales as short 
as 22 min. Some of the V404 optical flares lag >20 minutes 
behind the X-rays (Fig. 2a), and similar lags are also seen at 
mm-radio wavelengths. This may resemble the correlated X- 
ray/infrared/radio oscillations also referred to as “30-minutes" 
cycles of GRS 1915-1-105 (e.g. Fender & Pooley 1998; Mirabel 
et al. 1998; Rodriguez et al. 2008a). In GRS 1915-1-105, how¬ 
ever, these events are associated with hard X-ray dips preced¬ 
ing the flares and a clear softening at the X-ray peak, marking 


the disappearance of the Comptonization component (Rodriguez 
et al. 2008b). V404, on the contrary, remains hard even in the 
flaring states (Fig. 1), indicating a different mechanism respon¬ 
sible for the X-ray flaring (similar results were obtained from the 
1989 outburst, e.g. Zycki et al. 1999). One tempting possibility 
would be that the high-energy flares are direct boosted emission 
from a jet (blazar-like configuration). This would imply a jet axis 
not perpendicular to the orbital plane. Misaligned jets have been 
seen in GRO J1655-40 and V4641 Sgr (Maccarone 2002, and 
references therein). In the former, a rather modest Lorentz factor 
y~2.5 implies small relativistic boosting (Hjellming & Rupen 
1995). In the latter, y ranges from 10 up to 17, and the angle be¬ 
tween the jet axis and the orbital plane normal is as high as 50° 
(Maccarone 2002, and references therein). Significant Doppler 
boosting is expected in this case. 

Zycki et al. (1999) argued that the spectral and intensity vari¬ 
ability seen with Ginga in 1989 could be due to the evolution of a 
heavily absorbing medium. However, even with Nh >10^"^ cm“^ 
(Zycki et al. 1999), the activity above 20 keV is not affected by 
absorption, and hence the absorber alone cannot be responsible 
for the large variability we observe. The high energy flares could 
be due to the shock of the relativistic jets with the dense am¬ 
bient medium. Then optically thin synchrotron emission would 
be expected at X-ray energies, while our analysis favors thermal 
Comptonization. More simultaneous multi-wavelength observa¬ 
tions will help disentangling these different possibilities. 
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Table 1. List of the >6 Crab flares and their main properties. MJD 57193 is 2015 June 20*. 


Name 

StarC 

(MJD) 

Peak time 
(MJD) 

Stop" 

(MJD) 

CR 3-13 keV^ 
(cts/s) 

CR 2 O- 4 O keV* 
(cts/s) 

Properties" 

I 

57193.9217 

57193.9356 

57193.9402 

184 

1215 

Multiple 

II 

57193.9703 

57193.9981 

57194.0827 

181 

1234 

Multiple 

III 

57194.0827 

57194.1152 

57194.1428 

1055 

5209 

Complex 

IV 

57194.2232 

57194.3107 

57194.3938 

2010 

7040 

Isolated/Complex 

Va"' 

57194.6273 

57194.6399 

57194.6521 

473" 

3328 

Multiple/Complex, preceded by plateau 

Vb"' 

57194.6521 

57194.6579 

57194.6708 

852 

3999 

Multiple/Complex 

VI 

57194.6960 

57194.7346 

57194.7473 

129 

1974 

Isolated/Complex 

VII 

57194.9788 

57194.9996 

57195.0089 

459 

2200 

Multiple/Complex 

VIII 

57195.0089 

57195.0293 

57195.0501 

865 

4950 

Multiple/Complex 

IX 

57195.0582 

57195.0826 

57195.1095 

320 

2386 

Multiple/Complex 

X 

57195.2318 

57195.2503 

57195.2712-^ 

577« 

3472 

Multiple, preceded by succession of ~6 Crab peaks 

XI 

57195.4294 

57195.4388 

57195.4450 

857 

7036 

Multiple 

Xlla"' 

57195.4450 

57195.4573 

57195.4665 

401 

3525 

Multiple/Complex 

Xllb"' 

57195.4665 

57195.4723 

57195.4841 

1231 

6299 

Multiple, followed by plateau 

XIII 

57197.1373 

57197.1785 

57197.1924 

2076 

7081 

Multiple, preceded by plateau 

XIV 

57197.1924 

57197.2020 

57197.2067 

1240 

4368 

Multiple 

XV 

57197.2124 

57197.2228 

57197.2310 

210 

1793 

Multiple/Complex 

XVI 

57197.3450 

57197.3647 

57197.3705 

151 

1036 

Isolated 


“ Start (resp. stop) time of a flare is defined as the time 20^0 keV CR reaches 165 cts/s (1 Crab) during the increase (resp. 

decrease), or by the minimum reached before (resp. after) the increase (decrease) for multiple flares. The uncertainty on the 

times is +6x10“^ d. 

* Count rates at the peaks 

“Multiple" stands for series of well defined flares occurring in rapid repetition. “Complex" stands for flares showing various 
peaks. “Plateau" indicates a >1 Crab plateau. 

These peaks appear as single peaks in Fig. 1. They are in fact true multiples. 

The 3-13 keV peak time occurred about 200 s before the 20^0 keV one, indicating a potential hard lag. 
f Data gap at the end of the flare. The stop time is the last point before the gap. 

® The 3-13 keV peak time occurred about 200 s after the 20^0 keV one, indicating a potential soft lag. 
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Fig. 4. INTEGRAL LC?^of V404 over our ~4 day-long observations. All spectral domains considered for the LC extraction are shown here. The 
dashed line in the 20-40 keV panel represents the approximate level of Le;]^ we estimated. MJD 57193 is 2015 June 20*. 














































